CHAPTER 



Projectiles and 
Charged Particles 






In this chapter, I present two topics: the motion of projectiles subject to the forces of 
gravity and air resistance, and the motion of charged particles in uniform magnetic 
fields. Both problems lend themselves to solution using Newton’s laws in Cartesian 
coordinates, and both allow us to review and introduce some important mathematics. 
Above all, both are problems of great practical interest. 


2.1 Air Resistance 


Most introductory physics courses spend some time studying the motion of projectiles, 
but they almost always ignore air resistance. In many problems this is an excellent 
approximation; in others, air resistance is obviously important, and we need to know 
how to account for it. More generally, whether or not air resistance is significant, we 
need some way to estimate how important it really is. 

Let us begin by surveying some of the basic properties of the resistive force, or 
drag, f of the air, or other medium, through which an object is moving. (I shall 
generally speak of “air resistance” since air is the medium through which most 
projectiles move, but the same considerations apply to other gases and often to liquids 
as well.) The most obvious fact about air resistance, well known to anyone who rides 
a bicycle, is that it depends on the speed, v, of the object concerned. In addition, for 
many objects, the direction of the force due to motion through the air is opposite to 
the velocity v. For certain objects, such as a nonrotating sphere, this is exactly true, 
and for many it is a good approximation. You should, however, be aware that there 
are situations where it is certainly not true: The force of the air on an airplane wing 
has a large sideways component, called the lift, without which no airplanes could fly. 
Nevertheless, I shall assume that f and v point in opposite directions; that is, I shall 
consider only objects for which the sideways force is zero, or at least small enough 
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Figure 2.1 A projectile is subject to two forces, the force 
of gravity, w = mg, and the drag force of air resistance, 
f = -/(u)v. 


to be neglected. The situation is illustrated in Figure 2.1 and is summed up in the 
equation 


f = -f(v)y, (2.1) 

where v = v/ 1 v | denotes the unit vector in the direction of v, and f(v) is the magnitude 
of f. 

The function f(v ) that gives the magnitude of the air resistance varies with v in 
a complicated way, especially as the object’s speed approaches the speed of sound. 
However, at lower speeds it is often a good approximation to write 1 

f(v) = bv + cv 2 = / lin + / quad (2.2) 

where / lin and / quad stand for the linear and quadratic terms respectively, 

/lin and /quad • (2.3) 

The physical origins of these two terms are quite different: The linear term, / lin , arises 
from the viscous drag of the medium and is generally proportional to the viscosity of 
the medium and the linear size of the projectile (Problem 2.2). The quadratic term, 
/ quad , arises from the projectile’s having to accelerate the mass of air with which it is 
continually colliding; / quad is proportional to the density of the medium and the cross- 
sectional area of the projectile (Problem 2.4). In particular, for a spherical projectile 
(a cannonball, a baseball, or a drop of rain), the coefficients b and c in (2.2) have the 
form 


b = fiD and c = y D 2 (2.4) 

where D denotes the diameter of the sphere and the coefficients /S and y depend 
on the nature of the medium. For a spherical projectile in air at STP, they have the 
approximate values 


P = 1.6 x 10 4 N-s/m 2 


(2.5) 


1 Mathematically, Equation (2.2) is, in a sense, obvious. Any reasonable function is expected to 
have a Taylor series expansion, f — a + bv + cv 2 + • • •. For low enough v, the first three terms 
should give a good approximation, and, since / = 0 when v = 0 the constant term, a, has to be zero. 
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and 


y = 0.25 N-s 2 /m 4 . 


( 2 . 6 ) 


(For calculation of these two constants, see Problems 2.2 and 2.4.) You need to 
remember that these values are valid only for a sphere moving through air at STP. 
Nevertheless, they give at least a rough idea of the importance of the drag force even 
for nonspherical bodies moving through different gases at any normal temperatures 
and pressures. 

It often happens that we can neglect one of the terms in (2.2) compared to the other, 
and this simplifies the task of solving Newton’s second law. To decide whether this 
does happen in a given problem, and which term to neglect, we need to compare the 
sizes of the two terms: 



y D 
- v 

P 


f 1.6 x 10 3 



Dv 


(2.7) 


if we use the values (2.5) and (2.6) for a sphere in air. In a given problem, we have 
only to substitute the values of D and v into this equation to find out if one of the 
terms can be neglected, as the following example illustrates. 


example 2.1 A Baseball and Some Drops of Liquid 

Assess the relative importance of the linear and quadratic drags on a baseball 
of diameter D — 1 cm, traveling at a modest v = 5 m/s. Do the same for a drop 
of rain (D — 1 mm and v = 0.6 m/s) and for a tiny droplet of oil used in the 
Millikan oildrop experiment (D — 1.5 /rra and v = 5 x 10 -5 m/s). 

When we substitute the numbers for the baseball into (2.7) (remembering to 
convert the diameter to meters), we get 

- — % 600 [baseball], (2.8) | 

/lin ) 

| 

For this baseball, the linear term is clearly negligible and we need consider only 
the quadratic drag. If the ball is traveling faster, the ratio / quad // lin is even greater. 

At slower speeds the ratio is less dramatic, but even at 1 m/s the ratio is 100. In 
fact if v is small enough that the linear term is comparable to the quadratic, both 
terms are so small as to be negligible. Thus, for baseballs and similar objects, it 
is almost always safe to neglect / Un and take the drag force to be 

f = -ct 2 v. (2.9) [ 

For the raindrop, the numbers give 

— ~ 1 [raindrop]. (2.10) I 

/lin ] 

Thus for this raindrop the two terms are comparable and neither can be ne¬ 
glected — which makes solving for the motion more difficult. If the drop were 
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a lot larger or were traveling much faster, then the linear term would be negligi- 

i 

ble; and if the drop were much smaller or were traveling much slower, then the 
quadratic term would be negligible. But in general, with raindrops and similar 
objects, we are going to have to take both / lin and / ad into account. 

For the oildrop in the Millikan experiment the numbers give 

| 

| 10- 7 [Millikan oildrop], (2.11) 

/lin 

In this case, the quadratic term is totally negligible, and we can take 

*? 

f _ —bvx ■— —bx, (2.12) 


where the second, very compact form follows because, of course, vx = x. 


The moral of this example is clear: First, there are objects for which the drag force 
is dominantly linear, and the quadratic force can be neglected — notably, very small 
liquid drops in air, but also slightly larger objects in a very viscous fluid, such as a ball 
bearing moving through molasses. On the other hand, for most projectiles, such as golf 
balls, cannonballs, and even a human in free fall, the dominant drag force is quadratic, 
and we can neglect the linear term. This situation is a little unlucky because the linear 
problem is much easier to solve than the quadratic. In the following two sections, 
I shall discuss the linear case, precisely because it is the easier one. Nevertheless, it 
does have practical applications, and the mathematics used to solve it is widely used in 
many fields. In Section 2.4,1 shall take up the harder but more usual case of quadratic 
drag. 

To conclude this introductory section, I should mention the Reynolds number, an 
important parameter that features prominently in more advanced treatments of motion 
in fluids. As already mentioned, the linear drag / lin can be related to the viscosity of the 
fluid through which our projectile is moving, and the quadratic term / quad is similarly 
related to the inertia (and hence density) of the fluid. Thus one can relate the ratio 
/quad//lin to the fundamental parameters t], the viscosity, and q, the density, of the 
fluid (see Problem 2.3). The result is that the ratio / quad //ij n is of roughly the same 
order of magnitude as the dimensionless number R = Dvq/t), called the Reynolds 
number. Thus a compact and general way to summarize the foregoing discussion is 
to say that the quadratic drag / quad is dominant when the Reynolds number R is large, 
whereas the linear drag dominates when R is small. 


2.2 Linear Air Resistance 


Let us consider first a projectile for which the quadratic drag force is negligible, so 
that the force of air resistance is given by (2.12). We shall see directly that, because 
the drag force is linear in v, the equations of motion are very simple to solve. The two 
forces on the projectile are the weight w = mg and the drag force f = —bx, as shown 
in Figure 2.2. Thus the second law, mr = F, reads 


mr = mg — bx. 


(2.13) 
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Figure 2.2 The two forces on a projectile for which the 
force of air resistance is linear in the velocity, f = —b\. 


An interesting feature of this form is that, because neither of the forces depends on r, 
the equation of motion does not involve r itself (only the first and second derivatives 
of r). In fact, we can rewrite r as v, and (2.13) becomes 

mv = mg — b\, (2.14) 

a first-order differential equation for v. This simplification comes about because the 
forces depend only on v and not r. It means we have to solve only a first-order 
differential equation for v and then integrate v to find r. 

Perhaps the most important simplifying feature of linear drag is that the equation 
of motion separates into components especially easily. For instance, with x measured 
to the right and y vertically downward, (2.14) resolves into 

mv x = —bv x (2.15) 


and 


mVy=mg — bv y . (2.16) 

That is, we have two separate equations, one for v x and one for v y ; the equation for v x 
does not involve v y and vice versa. It is important to recognize that this happened only 
because the drag force was linear in v. For instance, if the drag force were quadratic. 



2 ^ 
CV V 


— cv\ — —c 



+ V 2 V, 

V 


(2.17) 


then in (2.14) we would have to replace the term — bv with (2.17). In place of the two 
equations (2.15) and (2.16), we would have 


mv x = —cJ v x + Vy v x 


mv y — mg — c./u 2 + u 2 Vy. 


y 


X 


y y 


(2.18) 


Here, each equation involves both of the variables v x and v y . These two coupled 
differential equations are much harder to solve than the uncoupled equations of the 
linear case. 

Because they are uncoupled, we can solve each equation for linear drag separately 
and then put the two solutions together. Further, each equation defines a problem that 
is interesting in its own right. Equation (2.15) is the equation of motion for an object 
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v 


f = -b\ 



Figure 2.3 A cart moves on a horizontal frictionless track 
in a medium that produces a linear drag force. 


(a cart with frictionless wheels, for instance) coasting horizontally in a medium that 
causes linear drag. Equation (2.16) describes an object (a tiny oil droplet for instance) 
that is falling vertically with linear air resistance. I shall solve these two separate 
problems in turn. 


Horizontal Motion with Linear Drag 

Consider an object such as the cart in Figure 2.3 coasting horizontally in a linearly 
resistive medium. I shall assume that at t = 0 the cart is at x =0 with velocity v x = v x0 . 
The only force on the cart is the drag f = — b\, thus the cart inevitably slows down. 
The rate of slowing is determined by (2.15), which has the general form 

v x = —kv x , (2.19) 

where k is my temporary abbreviation for k = b/m. This is a first-order differential 
equation for v x , whose general solution must contain exactly one arbitrary constant. 
The equation states that the derivative of v x is equal to —k times v x itself, and the only 
function with this property is the exponential function 

v x (t) = Ae~ kt (2.20) 

which satisfies (2.19) for any value of the constant A (Problems 1.24 and 1.25). Since 
this solution contains one arbitrary constant, it is the general solution of our first- 
order equation; that is, any solution must have this form. In our case, we know that 
v x (0) = v x0 , so that A = v x0 , and we conclude that 

v x (t) = v xo e~ kt = v xo e~ t/r , (2.21) 

where I have introduced the convenient parameter 

r = 1/Jfc = m/b [for linear drag]. (2.22) 

We see that our cart slows down exponentially, as shown in Figure 2.4(a). The 
parameter r has the dimensions of time (as you should check), and you can see from 

(2.21) that when t = t, the velocity is 1/e of its initial value; that is, t is the “1/e” 
time for the exponentially decreasing velocity. As t —> oo, the velocity approaches 
zero. 

To find the position as a function of time, we have only to integrate the velocity 

(2.21) . Integrations of this kind can be done using the definite or indefinite integral. 
The definite integral has the advantage that, it automatically takes care of the constant 
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Figure 2.4 (a) The velocity v x as a function of time, t, for a cart 

moving horizontally with a linear resistive force. As t —► oo, 14 - 
approaches zero exponentially, (b) The position x as a function of 
t for the same cart. As t -» 00 , x —> x^ = v xo x. 


of integration: Since v x = dx/dt, 

t 

v x (t') dt' = x(t) — x( 0 ). 

(Notice that I have named the “dummy” variable of integration t' to avoid confusion 
with the upper limit t .) Therefore 



x(f)=x(0) + / v xo e 

Jo 



= 0+ [-v xo re ,,/T ] o 

= x OQ (l- e ~ t/x ). (2.23) 

In the second line, I have used our assumption that x = 0 when t = 0. And in the last, 
I have introduced the parameter 


Xqq (2.24) 

which is the limit of x(t) as t —> 00 . We conclude that, as the cart slows down, its 
position approaches x^ asymptotically, as shown in Figure 2.4(b). 


Vertical Motion with Linear Drag 

Let us next consider a projectile that is subject to linear air resistance and is thrown 
vertically downward. The two forces on the projectile are gravity and air resistance, as 
shown in Figure 2.5. If we measure y vertically down, the only interesting component 
of the equation of motion is the y component, which reads 

miy — mg — bv y . (2.25) 

With the velocity downward (v y > 0), the retarding force is upward, while the force 
of gravity is downward. If v y is small, the force of gravity is more important than 
the drag force, and the falling object accelerates in its downward motion. This will 
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Figure 2.5 The forces on a projectile that is thrown ver¬ 
tically down, subject to linear air resistance. 


continue until the drag force balances the weight. The speed at which this balance 
occurs is easily found by setting (2.25) equal to zero, to give v y = mg/b or 


Vy = V 


ter 


where I have defined the terminal speed 



mg 

b 


[for linear drag]. 


(2.26) 


The terminal speed is the speed at which our projectile will eventually fall, if given 
the time to do so. Since it depends on m and b, it is different for different bodies. For 
example, if two objects have the same shape and size (b the same for both), the heavier 
object (m larger) will have the higher terminal speed, just as you would expect. Since 
u ter is inversely proportional to the coefficient b of air resistance, we can view v ter as 
an inverse measure of the importance of air resistance — the larger the air resistance, 
the smaller v ter , again just as you would expect. 


example 2.2 Terminal Speed of Small Liquid Drops 

•s 

j Find the terminal speed of a tiny oildrop in the Millikan oildrop experiment 

j (diameter D = 1.5 /xm and density q = 840 kg/m 3 ). Do the same for a small 

\ 

j drop of mist with diameter D = 0.2 mm. 

| 

j From Example 2.1 we know that the linear drag is dominant for these objects, 

I 

j so the terminal speed is given by (2.26). According to (2.4), b = /3D where 
j fi = 1.6 x 10~ 4 (in SI units). The mass of the drop is m = qtcD^/6. Thus (2.26) 
1 becomes 

I 

I 

j Q7tD Z g 

1 ^ter 

1 


6/J 


[for linear drag]. 


(2.27) 
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This interesting result shows that, for a given density, the terminal speed is 
proportional to D 2 . This implies that, once air resistance has become important, 
a large sphere will fall faster than a small sphere of the same density. 2 
Putting in the numbers, we find for the oildrop 


(840) XTr x (1.5 x 10~ 6 ) x (9.8) , , 1A _ 5 . 

v tPr = ----- i- --—- =6.1 x 10 5 m/s 

ter 6 x (1.6 x 10~ 4 ) 


[oildrop]. 


In the Millikan oildrop experiment, the oildrops fall exceedingly slowly, so their 
speed can be measured by simply watching them through a microscope. 
Putting in the numbers for the drop of mist, we find similarly that 


i> ter = 1.3 m/s [drop of mist]. 


(2.28) 


This speed is representative for a fine drizzle. For a larger raindrop, the terminal 
speed would be appreciably larger, but with a larger (and hence also faster) drop, 
the quadratic drag would need to be included in the calculation to get a reliable j 
value for t> ter . j 


So far, we have discussed the terminal speed of a projectile (moving vertically), 
but we must now discuss how the projectile approaches that speed. This is determined 
by the equation of motion (2.25) which we can rewrite as 


mv y = -b(v y - u ter ). (2.29) 

(Remember that v [er — mg/b.) This differential equation can be solved in several ways. 
(For one alternative see Problem 2.9.) Perhaps the simplest is to note that it is almost 
the same as Equation (2.15) for the horizontal motion, except that on the right we now 
have (v y — u ter ) instead of v x . The solution for the horizontal case was the exponential 
function (2.20). The trick to solving our new vertical equation (2.29) is to introduce 
the new variable u = (v y — u ter ), which satisfies mu = —bu (because u ter is constant). 
Since this is exactly the same as Equation (2.15) for the horizontal motion, the solution 
for u is the same exponential, u = Ae~^ x . [Remember that the constant k in (2.20) 
became k = 1/r.] Therefore, 

v y ~ v ter = Ae~ t/z . 

When t = 0, v y = v y0 , so A = v y0 — u ter and our final solution for v y as a function of 
t is 

Vy(t) = V ter + (v y0 - v ter )e~ t/T (2.30) 

= v yo e~ t/r + v ter (l - e~ t/T ) . (2.31) 


2 We are here assuming that the drag force is linear, but the same qualitative conclusion follows 
for a quadratic drag force. (Problem 2.24.) 
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Figure 2.6 When an object is dropped in a medium with 
linear resistance, v y approaches its terminal value v ter as 
shown. 


This second expression gives v y (t ) as the sum of two terms: The first is equal to v y0 
when t = 0, but fades away to zero as t increases; the second is equal to zero when 
t = 0, but approaches u ter as t —> oo. In particular, as t —> oo, 

Vyit) -* u ter (2.32) 

just as we anticipated. 

Let us examine the result (2.31) in a little more detail for the case that v y0 = 0; 
that is, the projectile is dropped from rest. In this case (2.31) reads 

v y (t) = v ter (l - e~ t/x ) . (2.33) 

This result is plotted in Figure 2.6, where we see that v y starts out from 0 and 
approaches the terminal speed, v y v teT , asymptotically as t -» oo. The significance 
of the time r for a falling body is easily read off from (2.33). When t = r, we see that 

v y = %r(l - e~ l ) = 0.63u ter . 

That is, in a time r, the object reaches 63% of the terminal speed. Similar calculations 
give the following results; 

time percent 

* °f %er 

0 0 

r 63% 

2t 86% 

3t 95% 

Of course, the object’s speed never actually reaches u ter , but r is a good measure of 
how fast the speed approaches u ter . In particular, when t — 3r the speed is 95% of u ter , 
and for many purposes we can say that after a time 3 t the speed is essentially equal 
to Uj er . 
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example 2.3 Characteristic Time for Two Liquid Drops 

Find the characteristic times, r, for the oildrop and drop of mist in Example 2.2. 

The characteristic time r was defined in (2.22) as r = m /b, and t> ter was j 
defined in (2.26) as n ter = mg lb. Thus we have the useful relation 

v ter = gx. (2.34) j 

Notice that this relation lets us interpret u ter as the speed a falling object would 
acquire in a time x, if it had a constant acceleration equal to g. Also note that, 
like u ter , the time r is an inverse indicator of the importance of air resistance: 
When the coefficient b of air resistance is small, both v ter and r are large; when 
b is large, both u ter and r are small. 

For our present purposes, the importance of (2.34) is that, since we have 
already found the terminal velocities of the two drops, we can immediately find 
the values of r. For the Millikan oildrop, we found that u ter = 6.1 x 10 -5 m/s, 
therefore j 

x = ^SL = 6 - X - 1Q 5 - 6.2 x 1CT 6 s [oildrop], 

8 9.8 

After falling for just 20 microseconds, this oildrop will have acquired 95% of 
its terminal speed. For almost every purpose, the oildrop always travels at its 
terminal speed. ) 

For the drop of mist of Example 2.2, the terminal speed was v ter — 1.3 m/s 
and so x = u ter /g ~ 0.13 s. After about 0.4 s, the drop will have acquired 95% 
of its terminal speed. 


Whether or not our falling object starts from rest, we can find its position y as a 
function of time by integrating the known form (2.30) of v y , 

0 %er "F (ny 0 Uj er )c 

Assuming that the projectile’s initial position is y = 0, it immediately follows that 



= v teT t + ( v y0 - u ter )r (1-e t/T ). (2.35) 

This equation for y(t) can now be combined with Equation (2.23) for x(t) to give 
us the orbit of any projectile, moving both horizontally and vertically, in a linear 
medium. 
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2.3 Trajectory and Range in a Linear Medium 


We saw at the begining of the last section that the equation of motion for a projectile 
moving in any direction resolves into two separate equations, one for the horizontal 
and one for the vertical motion [Equations (2.15) and (2.16)]. We have solved each 
of these separate equations in (2.23) and (2.35), and we can now put these solutions 
together to give the trajectory of an arbitrary projectile moving in any direction. In 
this discussion it is marginally more convenient to measure y vertically upward, in 
which case we must reverse the sign of u ter . (Make sure you understand this point.) 
Thus the two equations of the orbit become 

x(t) = v X0 x (l — e~^ z ) 1 

y(t) = (v y0 + u ter )r (l - e t,r ) - v teT t. 

You can eliminate t from these two equations by solving the first for t and then 
substituting into the second. (See Problem 2.17.) The result is the equation for the 
trajectory: 


v yo + ^ter , , f i X \ /O 'in'i 

y =-x + i> ter r In I 1-1 . (2.37) 

Lto V v xo r J 

This equation is probably too complicated to be especially illuminating, but I have 
plotted it as the solid curve in Figure 2.7, with the help of which you can understand 
some of the features of (2.37). For example, if you look at the second term on the right 
of (2.37), you will see that as x -> v xo r the argument of the log function approaches 
zero; therefore, the log term and hence y both approach —oo. That is, the trajectory 
has a vertical asymptote at x = v xo x, as you can see in the picture. I leave it as an 
exercise (Problem 2.19) for you to check that if air resistance is switched off (u ter and 
x both approach infinity), the trajectory defined by (2.37) does indeed approach the 
dashed trajectory corresponding to zero air resistance. 


Horizontal Range 

A standard (and quite interesting) problem in elementary physics courses is to show 
that the horizontal range R of a projectile (subject to no air resistance of course) is 

2v v 

i? vac = — -° - v ° [no air resistance] (2.38) 

8 

where R vac stands for the range in a vacuum. Let us see how this result is modified by 
air resistance. 

The range R is the value of x when y as given by (2.37) is zero. Thus R is the 
solution of the equation 


— R + twr In (l -— ^ 

V v xo x) 



v yo + V \ 


(2.39) 
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► 


Figure 2.7 The trajectory of a projectile subject to a linear drag 
force (solid curve) and the corresponding trajectory in a vacuum 
(dashed curve). At first the two curves are very similar, but as t 
increases, air resistance slows the projectile and pulls its trajec¬ 
tory down, with a vertical asymptote at x = v xo z. The horizontal 
range of the projectile is labeled R, and the corresponding range 
in vacuum /? vac . 


This is a transcendental equation and cannot be solved analytically, that is, in terms of 
well known, elementary functions such as logs, or sines and cosines. For a given choice 
of parameters, it can be solved numerically with a computer (Problem 2.22), but this 
approach usually gives one little sense of how the solution depends on the parameters. 
Often a good alternative is to find some approximation that allows an approximate 
analytic solution. (Before the advent of computers, this was often the only way to find 
out what happens.) In the present case, it is often clear that the effects of air resistance 
should be small. This means that both u ter and r are large and the second term in the 
argument of the log function is small (since it has r in its denominator). This suggests 
that we expand the log in a Taylor series (see Problem 2.18): 

ln(l f) = — H - • • • (2.40) 

We can use this expansion for the log term in (2.39), and, provided r is large enough, 
we can surely neglect the terms beyond e 3 . This gives the equation 


VyO + U[ er 


V 


XO 


^ ^ter^ 


2L + 2 (2L 

V X o r 2 \ v xo T 





(2.41) 


This equation can be quickly tidied up. First, the second term in the first bracket 
cancels the first term in the second. Next, every term contains a factor of R. This 
implies that one solution is R = 0, which is correct — the height y is zero when x = 0. 
Nevertheless, this is not the solution we are interested in, and we can divide out the 
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common factor of R. A little rearrangement (and replacement of u ter /r by g) lets us 
rewite the equation as 


2v xo v vo 2 ^ 

R = - - f . 2 y° - — R 2 . (2.42) 

8 3u xo r 

This may seem a perverse way to write a quadratic equation for R, but it leads us 
quickly to the desired approximate solution. The point is that the second term on the 
right is very small. (In the numerator R is certainly no more than 7? vac and we are 
assuming that r in the denominator is very large.) Therefore, as a first approximation 
we get 


R 




8 



(2.43) 


This is just what we expected: For low air resistance, the range is close to /? vac . But 
with the help of (2.42) we can now get a second, better approximation. The last term of 
(2.42) is the required correction to i? vac ; because it is already small, we would certainly 
be satisfied with an approximate value for this correction. Thus, in evaluating the last 
term of (2.42), we can replace R with the approximate value R ~ i? vac , and we find 
as our second approximation [remember that the first term in (2.42) is just /? vae ] 


R ~ tfvac - --C^vac ) 2 

3 v xo r 

= ^vac(l-^—)• (2-44) 

V 3 Uj er ) 

(To get the second line, I replaced the second R vac in the previous line by 2v xo v yo /g 
and r g by i> ter .) Notice that the correction for air resistance always makes R smaller 
than /? vac , as one would expect. Notice also that the correction depends only on the 
ratio v y0 /v tN . More generally, it is easy to see (Problem 2.32) that the importance of 
air resistance is indicated by the ratio v/v ter of the projectile’s speed to the terminal 
speed. If v/v ter <<C 1 throughout the flight, the effect of air resistance is very small; 
if v/v ter is around 1 or more, air resistance is almost certainly important [and the 
approximation (2.44) is certainly no good]. 


example 2.4 Range of Small Metal Pellets 

I flick a tiny metal pellet with diameter d = 0.2 mm and v = 1 m/s at 45°. Find 
its horizontal range assuming the pellet is gold (density q 16g/cm 3 ). What if 
it is aluminum (density q ~ 2.7 g/cm 3 )? 

In the absence of air resistance, both pellets would have the same range, 

D 2v xo v 0 

R wac = -= 10.2 cm. 


g 
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For gold, Equation (2.27) gives (as you can check) v ler ~ 21 m/s. Thus the 
correction term in (2.44) is 


4 v 4 0.71 

= - x - 

3 21 


3 v 


0.05. 


ter 


That is, air resistance reduces the range by 5% to about 9.7 cm. The density 
of aluminum is about 1/6 times that of gold. Therefore the terminal speed is 
one sixth as big, and the correction for aluminum is 6 times greater or about 
30%, giving a range of about 7 cm. For the gold pellet the correction for air 
resistance is quite small and could perhaps be neglected; for the aluminum pellet, 
the correction is still small, but is certainly not negligible. 


2.4 Quadratic Air Resistance 


In the last two sections we have developed a rather complete theory of projectiles 
subject to a linear drag force, f = —by. While we can find examples of projectiles for 
which the drag is linear (notably very small objects, such as the Millikan oildrop), for 
most of the more obvious examples of projectiles (baseballs, footballs, cannonballs, 
and the like) it is a far better approximation to say that the drag is pure quadratic, 
f = —cv 2 \. We must, therefore, develop a corresponding theory for a quadratic drag 
force. On the face of it, the two theories are not so very different. In either case we 
have to solve the differential equation 

m\ = mg + f, (2.45) 

and in both cases this is a first-order differential equation for the velocity v, with f 
depending in a relatively simple way on v. There is, however, an important difference. 
In the linear case (f = —by). Equation (2.45) is a linear differential equation, inas¬ 
much as the terms that involve v are all linear in v or its derivatives. In the quadratic 
case, Equation (2.45) is, of course, nonlinear. And it turns out that the mathemati¬ 
cal theory of nonlinear differential equations is significantly more complicated than 
the linear theory. As a practical matter, we shall find that for the case of a general 
projectile, moving in both the x and y directions. Equation (2.45) cannot be solved 
in terms of elementary functions when the drag is quadratic. More generally, we 
shall see in Chapter 12 that for more complicated systems, nonlinearity can lead to 
the astonishing phenomenon of chaos, although this does not happen in the present 
case. 

In this section, I shall start with the same two special cases discussed in Section 2.2, 
a body that is constrained to move horizontally, such as a railroad car on a horizontal 
track, and a body that moves vertically, such as a stone dropped from a window (both 
now with quadratic drag forces). We shall find that in these two especially simple cases 
the differential equation (2.45) can be solved by elementary means, and the solutions 
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introduce some important techniques and interesting results. I shall then discuss briefly 
the general case (motion in both the horizontal and vertical directions), which can be 
solved only numerically. 


Horizontal Motion with Quadratic Drag 

Let us consider a body moving horizontally (in the positive x direction), subject to 
a quadratic drag and no other forces. For example, you could imagine a cycle racer, 
who has crossed the finishing line and is coasting to a stop under the influence of air 
resistance. To the extent that the cycle is well lubricated and tires well inflated, we 
can ignore ordinary friction, 3 and, except at very low speeds, air resistance is purely 
quadratic. The x component of the equation of motion is therefore (I’ll abbreviate v x 
to v ) 


dv 
m — 
dt 



(2.46) 


If we divide by v 2 and multiply by dt, we get an equation in which only the variable 
v appears on the left and only t on the right: 4 

= -cdt. (2.47) 

This trick — of rearranging a differential equation so that only one variable appears 
on the left and only the other on the right — is called separation of variables. When 
it is possible, separation of variables is often the simplest way to solve a first-order 
differential equation, since the solution can be found by simple integration of both 
sides. 

Integrating Equation (2.47) we find 



where v 0 is the initial velocity at t = 0. Notice that I have written both sides as definite 
integrals, with the appropriate limits, so that I shan’t have to worry about any constants 
of integration. I have also renamed the variables of integration as v' and /' to avoid 


3 As I shall discuss shortly, when the cyclist slows down to a stop, air resistance becomes smaller, 
and eventually friction becomes the dominant force. Nevertheless, at speeds around 10 mph or more, 
it is a fair approximation to ignore everything but the quadratic air resistance. 

4 In passing from (2.46) to (2.47), I have treated the derivative dv/dt as if it were the quotient of 
two separate numbers, dv and dt. As you are certainly aware this cavalier proceeding is not strictly 
correct. Nevertheless, it can be justified in two ways. First, in the theory of differentials , it is in fact 
true that dv and dt are defined as separate numbers (differentials), such that their quotient is the 
derivative dv/dt. Fortunately, it is quite unnecessary to know about this theory. As physicists we 
know that dv/dt is the limit of Av/ At, as both Av and At become small, and I shall take the view 
that dv is just shorthand for Av (and likewise dt for At), with the understanding that it has been 
taken small enough that the quotient dv / dt is within my desired accuracy of the true derivative. With 
this understanding, (2.47), with dv on one side and dt on the other, makes perfectly good sense. 
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confusion with the upper limits v and t. Both of these integrals are easily evaluated, 
and we find 



(2.48) 


or, solving for v, 


v(t) = - 2 -=-(2.49) 

1 + cv 0 t jm 

where I have introduced the abbreviation r for the combination of constants 

771 

r = — [for quadratic drag]. (2.50) 

cv 0 

As you can easily check, r is a time, with the significance that when t — r the velocity 
is v = v 0 /2. Notice that this parameter r is different from the x introduced in (2.22) for 
motion subject to linear air resistance; nevertheless, both parameters have the same 
general significance as indicators of the time for air resistance to slow the motion 
appreciably. 

To find the bicycle’s position jc, we have only to integrate v to give (as you should 
check) 


x(t)=jc 0 + f v(t')dt' 

Jo 

= v 0 r In (1 + t/x ), (2.51) 

if we take the initial position x 0 to be zero. Figure 2.8 shows our results for v and 
x as functions of t. It is interesting to compare these graphs with the corresponding 
graphs of Figure 2.4 for a body coasting horizontally but subject to a linear resistance. 
Superficially, the two graphs for the velocity look similar. In particular, both go to 
zero as t —* oo. But in the linear case v goes to zero exponentially, whereas in the 
quadratic case it does so only very slowly, like 1 ft. This difference in the behavior 
of v manifests itself quite dramatically in the behavior of x. In the linear case, we 



Figure 2.8 The motion of a body, such as a bicycle, coasting 
horizontally and subject to a quadratic air resistance, (a) The 
velocity is given by (2.49) and goes to zero like \/t as t — » oo. 
(b) The position is given by (2.51) and goes to infinity as t — > oo. 
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saw that x approaches a finite limit as t oo, but it is clear from (2.51) that in the 
quadratic case x increases without limit as t -* oo. 

The striking difference in the behavior of x for quadratic and linear drags is easy 
to understand qualitatively. In the quadratic case, the drag is proportional to v 2 . Thus 
as v gets small, the drag gets very small — so small that it fails to bring the bicyle 
to rest at any finite value of x. This unexpected behavior serves to highlight that a 
drag force that is proportional to v 2 at all speeds is unrealistic. Although the linear 
drag and ordinary friction are very small, nevertheless as v —> 0 they must eventually 
become more important than the v 2 term and cannot be ignored. In particular, one or 
another of these two terms (friction in the case of a bicycle) ensures that no real body 
can coast on to infinity! 


Vertical Motion with Quadratic Drag 

The case that an object moves vertically with a quadratic drag force can be solved in 
much the same way as the horizontal case. Consider a baseball that is dropped from a 
window in a high tower. If we measure the coordinate y vertically down, the equation 
of motion is (I’ll abbreviate v y to v now) 

mv = mg — cir. (2.52) 

Before we solve this equation, let us consider the ball’s terminal speed, the speed at 
which the two terms on the right of (2.52) just balance. Evidently this must satisfy 
cv 2 — mg, whose solution is 


(2.53) 

For any given object (given m, g, and c), this lets us calculate the terminal speed. For 
example, for a baseball it gives (as we shall see in a moment) u ter ~ 35 m/s, or nearly 
80 miles per hour. 

We can tidy the equation of motion (2.52) a little by using (2.53) to replace c by 
mg Jv 2 and canceling the factors of m: 



v = g 




(2.54) 


This can be solved by separation of variables, just as in the case of horizontal 
motion: First we can rewrite it as 



(2.55) 


This is the desired separated form (only v on the left and only t on the right) and 
we can simply integrate both sides. 5 Assuming the ball starts from rest, the limits of 


5 Notice that in fact any one-dimensional problem where the net force depends only on the 
velocity can be solved by separation of variables, since the equation mv = F(v) can always be 
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integration are 0 and v on the left and 0 and t on the right, and we find (as you should 
verify — Problem 2.35) 


arctanh (——) = t (2.56) 

g V / 

where “arctanh” denotes the inverse hyperbolic tangent. This particular integral can be 
evaluated alternatively in terms of the natural log function (Problem 2.37). However, 
the hyperbolic functions, sinh, cosh, and tanh, and their inverses arcsinh, arccosh, and 
arctanh, come up so often in all branches of physics that you really should learn to use 
them. If you have not had much exposure to them, you might want to look at Problems 
2.33 and 2.34, and study graphs of these functions. 

Equation (2.56) can be solved for v to give 


v — u ter tanh 


(-)' 

V ^ter / 


(2.57) 


To find the position y, we just integrate v to give 


y= (Uter)_ ln 
g 


cosh 


gt 


V 


ter/ J 


(2.58) 


While both of these two formulas can be cleaned up a little (see Problem 2.35), they 
are already sufficient to work the following example. 


example 2.5 A Baseball Dropped from a High Tower 

Find the terminal speed of a baseball (mass m = 0.15 kg and diameter D = 7 
cm). Make plots of its velocity and position for the first six seconds after it is 
dropped from a tall tower. 

The terminal speed is given by (2.53), with the coefficient of air resistance c 
given by (2.4) as c = y D 2 where y — 0.25 N-s 2 /m 4 . Therefore 

„ ter = = J (0-lSkg) X (9.8m/s 2 ) 35m/s (2 59 j 

V yD 2 V (0.25N-s 2 /m 4 ) x (0.07 m) 2 

or nearly 80 miles per hour. It is interesting to note that fast baseball pitchers 
can pitch a ball considerably faster than u ter . Under these conditions, the drag 
force is actually greater than the ball’s weight! 

The plots of v and y can be made by hand, but are, of course, much easier with 
the help of computer software such as Mathcad or Mathematica that can make the 
plots for you. Whatever method we choose, the results are as shown in Figure 2.9, 
where the solid curves show the actual velocity and position while the dashed 
curves are the corresponding values in a vacuum. The actual velocity levels out, I 


written as m dv/F(v) = dt. Of course there is no assurance that this can be integrated analytically 
if F(y) is too complicated, but it does guarantee a straightforward numerical solution at worst. See 
Problem 2.7. 
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(a) (b) 

Figure 2.9 The motion of a baseball dropped from the top of a high tower (solid 
curves). The corresponding motion in a vacuum is shown with long dashes, 
(a) The actual velocity approaches the ball’s terminal velocity r tcr = 35 m/s as 
t —> oo. (b) The graph of position against time falls further and further behind 
the corresponding vacuum graph. When t — 6 s, the baseball has dropped about 
130 meters; in a vacuum, it would have dropped about 180 meters. 


approaching the terminal value u ter = 35 m/s as t —> oo, whereas the velocity 
in a vacuum would increase without limit. Initially, the position increases just 
as it would in a vacuum (that is, y = \gt 2 ), but falls behind as v increases and 
j the air resistance becomes more important. Eventually, y approaches a straight 
line of the form y = v teT t+ const. (See Problem 2.35.) 


Quadratic Drag with Horizontal and Vertical Motion 


The equation of motion for a projectile subject to quadratic drag. 


mr = mg — cirv 
= mg — cvx, 


(2.60) 


resolves into its horizontal and vertical components (with y measured vertically 
upward) to give 


mv x 

mv y 


-CjV 2 + V 2 Vy 

V X ' y X 


m g ~ CjV 2 + V 2 Vy. 


(2.61) 


These are two differential equations for the two unknown functions v x (t) and v y (t), 
but each equation involves both v x and v y . In particular, neither equation is the same 
as for an object that moves only in the x direction or only in the y direction. This 
means that we cannot solve these two equations by simply pasting together our two 
separate solutions for horizontal and vertical motion. Worse still, it turns out that 
the two equations (2.61) cannot be solved analytically at all. The only way to solve 
them is numerically, which we can only do for specified numerical initial conditions 
(that is, specified values of the initial position and velocity). This means that we 
cannot find the general solution; all we can do numerically is to find the particular 
solution corresponding to any chosen initial conditions. Before I discuss some general 
properties of the solutions of (2.61), let us work out one such numerical solution. 
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example 2.6 Trajectory of a Baseball 

The baseball of Example 2.5 is now thrown with velocity 30 m/s (about 70 mi/h) j 
at 50° above the horizontal from a high cliff. Find its trajectory for the first eight j 

seconds of flight and compare with the corresponding trajectory in a vacuum. If j 

| 

the same baseball was thrown with the same initial velocity on horizontal ground 
how far would it travel before landing? That is, what is its horizontal range? 

We have to solve the two coupled differential equations (2.61) with the initial j 
conditions 


v x0 = v 0 cos 9 = 19.3 m/s and 


v yo = v o s i n $ = 23.0 m/s 


and x 0 = y Q = 0 (if we put the origin at the point from which the ball is thrown). 
This can be done with systems such as Mathematica, Matlab, or Maple, or with 
programming languages such as “C” or Fortran. Figure 2.10 shows the resulting 
trajectory, found using the function “NDSolve” in Mathematica. 

Several features of Figure 2.10 deserve comment. Obviously the effect of 
air resistance is to lower the trajectory, as compared to the vacuum trajectory 
(shown dashed). For example, we see that in a vacuum the high point of the 
trajectory occurs at f ^ 2.3 s and is about 27 m above the starting point; with air 
resistance, the high point comes just before t = 2.0 s and is at about 21 m. In 
a vacuum, the ball would continue to move indefinitely in the x direction. The 




Figure 2.10 Trajectory of a baseball thrown off a cliff and subject 
to quadratic air resistance (solid curve). The initial velocity is 30 
m/s at 50° above the horizontal; the terminal speed is 35 m/s. The 
dashed curve shows the corresponding trajectory in a vacuum. 
The dots show the ball’s position at one-second intervals. Air 
resistance slows the horizontal motion, so that the ball approaches 
a vertical asymptote just beyond x = 100 meters. 
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effect of air resistance is to slow the horizontal motion so that x never moves to 
the right of a vertical asymptote near x = 100 m. 

I 

The horizontal range of the baseball is easily read off the figure as the value 
of x when y returns to zero. We see that ~ 59 m, as opposed to the range in 
vacuum, R vac ~ 90 m. The effect of air resistance is quite large in this example, 
as we might have anticipated: The ball was thrown with a speed only a little 
j less than the terminal speed (30 vs 35 m/s), and this means that the force of air 
j resistance is only a little less than that of gravity. This being the case, we should 
expect air resistance to change the trajectory appreciably. 


This example illustrates several of the general features of projectile motion with a 
quadratic drag force. Although we cannot solve analytically the equations of motion 
(2.61) for this problem, we can use the equations to prove various general properties 
of the trajectory. For example, we noticed that the baseball reached a lower maximum 
height, and did so sooner, than it would have in a vacuum. It is easy to prove that this 
will always be the case: As long as the projectile is moving upward (v y > 0), the force 
of air resistance has a downward y component. Thus the downward acceleration is 
greater than g (its value in vacuum). Therefore a graph of v y against t slopes down from 
v y0 more quickly than it would in vacuum, as shown in Figure 2.11. This guarantees 
that v y reaches zero sooner than it would in vacuum, and that the ball travels less 
distance (in the y direction) before reaching the high point. That is, the ball’s high 
point occurs sooner, and is lower, than it would be in a vacuum. 



Figure 2.11 Graph of v y against t for a projectile that is thrown 
upward ( v yo > 0) and is subject to a quadratic resistance (solid 
curve). The dashed line (slope — — g) is the corresponding graph 
when there is no air resistance. The projectile moves upward 
until it reaches its maximum height when t> = 0. During this 
time, the drag force is downward and the downward acceleration 
is always greater than g. Therefore, the curve slopes more steeply 
than the dashed line, and the projectile reaches its high point 
sooner than it would in a vacuum. Since the area under the curve 
is less than that under the dashed line, the projectile’s maximum 
height is less than it would be in a vacuum. 
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I claimed that the baseball of Example 2.6 approaches a vertical asymptote as 
t oo, and we can now prove that this is always the case. First, it is easy to 
convince yourself that once the ball starts moving downward, it continues to accelerate 
downward, with v y approaching — u ter as t -» oo. At the same time v x continues to 
decrease and approaches zero. Thus the square root in both of the equations (2.61) 
approaches i> ter . In particular, when t is large, the equation for v x can be approximated 
by 



^ter 

m 


v 



say. The solution of this equation is, of course, an exponential function, v x = Ae~ kt , 
and we see that v x approaches zero very rapidly (exponentially) as t —► oo. This 
guarantees that x, which is the integral of v x . 


x(t) — f v x (t')dt', 

Jo 


approaches a finite limit as t —> oo, and the trajectory has a finite vertical asymptote 
as claimed. 


2.5 Motion of a Charge in a Uniform Magnetic Field 


Another interesting application of Newton’s laws, and (like projectile motion) an 
application that lets me introduce some important mathematical methods, is the 
motion of a charged particle in a magnetic field. I shall consider here a particle of 
charge q (which I shall usually take to be positive), moving in a uniform magnetic 
field B that points in the z direction as shown in Figure 2.12. The net force on the 
particle is just the magnetic force 


F = q\ x B, 


(2.62) 



Figure 2.12 A charged particle moving in a uniform mag¬ 
netic field that points in the z direction. 
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so the equation of motion can be written as 

mv = q\ x B. (2.63) 

[As with projectiles, the force depends only on the velocity (not the position), so the 
second law reduces to a first-order differential equation for v.] 

As is so often the case, the simplest way to solve the equation of motion is to 
resolve it into components. The components of v and B are 

V = (V x , Vy, V z ) 

and 

B = (0,0, B ), 

from which we can read off the components of v x B: 

v x B = ( v y B, —v x B, 0). 

Thus the three components of (2.63) are 

mv x = q B v y 
mv y = —qBv x 
mv z = 0. 

The last of these says simply that v z , the component of the 
direction of B, is constant: 


(2.64) 

(2.65) 

( 2 . 66 ) 

particle’s velocity in the 


v z = const, 

a result we could have anticipated since the magnetic force is always perpendicular to 
B. Because v z is constant, we shall focus most of our attention on v x and v y . In fact, 
we can even think of them as comprising a two-dimensional vector ( v x , v y ), which is 
just the projection of v onto the xy plane and can be called the transverse velocity , 


(v x , v y ) = transverse velocity. 

To simplify the equations (2.64) and (2.65) for v x and v y , I shall define the 
parameter 


co = 


qB 

m 


(2.67) 


which has the dimensions of inverse time and is called the cyclotron frequency. With 
this notation. Equations (2.64) and (2.65) become 


V X = (DVy 

Vy (OV X . 


( 2 . 68 ) 


These two coupled differential equations can be solved in a host of different ways. 
I would like to describe one that makes use of complex numbers. Though perhaps 
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real part 


Figure 2.13 The complex number rj = v x + iv y is repre¬ 
sented as a point in the complex plane. The arrow pointing 
from O to r) is literally a picture of the transverse velocity 
vector (v x , v y ). 


not the easiest solution, this method has surprisingly wide application in many ar¬ 
eas of physics. (For an alternative solution that avoids complex numbers, see Prob¬ 
lem 2.54.) 

The two variables v x and v y are, of course, real numbers. However, there is nothing 
to prevent us from defining a complex number 

r] = v x + iv y , (2.69) 

where i denotes the square root of — 1 (called j by most engineers), i — V—1 (and 
rj is the Greek letter eta). If we draw the complex number rj in the complex plane, 
or Argand diagram, then its two components are v x and v y as shown in Figure 2.13; 
in other words, the representation of r] in the complex plane is a picture of the two- 
dimensional transverse velocity ( v x , v y ). 

The advantages of introducing the complex number rj appear when we evaluate its 
derivative. Using (2.68), we find that 

rj — v x 4- iv y = cov y — ioov x — —ico(v x + iv y ) 


or 


rj = —icoij. (2.70) 

We see that the two coupled equations for v x and v have become a single equation for 
the complex number rj. Furthermore, it is an equation of the now familiar form u = ku, 
whose solution we know to be the exponential u = Ae kt . Thus we can immediately 
write down the solution for rj : 


rj = Ae~ im (2.71) 

Before we discuss the significance of this solution, I would like to review a few 
properties of complex exponentials in the next section. If you are very familiar with 
these ideas, by all means skip this material. 
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2.6 Complex Exponentials 


While you are certainly familiar with the exponential function e x for a real variable 
x, you may not be so at home with e z when z is complex. 6 For the real case there are 
several possible definitions of e x (for instance, as the function that is equal to its own 
derivative). The definition that extends most easily to the complex case is the Taylor 
series (see Problem 2.18) 


2 3 

e z =l + z + - + - + -- -. (2.72) 

2! 3! 

For any value of z, real or complex, large or small, this series converges to give a well- 
defined value for e z . By differentiating it, you can easily convince yourself that it has 
the expected property that it equals its own derivative. And one can show (not always 
so easily) that it has all the other familiar properties of the exponential function — for 
instance, that e z e w = e iz+w> . (See Problems 2.50 and 2.51.) In particular, the function 
Ae kz (with A and k any constants, real or complex) has the property that 

j- = k (a/ z ) . (2.73) 

Since it satisfies this same equation whatever the value of A, it is, in fact, the 
general solution of the first-order equation df/dz = kf. At the end of the last section, 
I introduced the complex number r](t) and showed that it satisfied the equation 
r) = —icor). We are now justified in saying that this guarantees that r] must be the 
exponential function anticipated in (2.71). 

We shall be particularly concerned with the exponential of a pure imaginary 
number, that is, e l ° where # is a real number. The Taylor series (2.72) for this function 
reads 


w , (i#) 2 (id) 3 (i#) 4 

10 = 1 + W + - + -—- + -— 3 _ + 


2! 


3! 


4! 


(2.74) 


Noting that i 2 — — 1, i 3 = —i, and so on, you can see that all of the even powers 
in this series are real, while all of the odd powers are pure imaginary. Regrouping 
accordingly, we can rewrite (2.74) to read 


e id = 


— 1 

<£> 

fO 

Q2> 

_1 


r # 3 l 

1-+ — + ••• 

+ i 

9 — — 

2! 4! J 

i 

LO 

1_ 


(2.75) 


The series in the first brackets is the Taylor series for cos#, and that in the second 
brackets is sin # (Problem 2.18). Thus we have proved the important relation: 


e' e — cos# + i sin#. 


(2.76) 


6 For a review of some elementary properties of complex numbers, see Problems 2.45 to 2.49. 
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Figure 2.14 (a) Euler’s formula, (2.76), implies that the complex num¬ 

ber e‘ e lies on the unit circle (the circle of radius 1, centered on the origin 
O) with polar angle 0. (b) The complex constant A = ae lS lies on a cir¬ 
cle of radius a with polar angle 8. The function rj(t) — Ae~ lcot lies on the 
same circle but with polar angle (S — cot) and moves clockwise around 
the circle as t advances. 


This result, known as Euler’s formula, is illustrated in Figure 2.14(a). Note especially 
that the complex number e' e has polar angle 6 and, since cos 2 9 + sin 2 # = 1, the 
magnitude of e‘° is 1; that is, e l ° lies on the unit circle, the circle with radius 1 centered 
at O. 

Our main concern is with a complex number of the form rj = Ae~ lC0t . The co¬ 
efficient A is a fixed complex number, which can be expressed as A = ae lS , where 
a — | A | is the magnitude, and <5 is the polar angle of A, as shown in Figure 2.14(b). 
(See Problem 2.45.) The number rj can therefore be written as 

r, = Ae~ ,0Jt = ae iS e- ia)t = ae i(S ~ wt) . (2.77) 


Thus r] has the same magnitude as A (namely a), but has polar angle equal to (8 — cot), 
as shown in Figure 2.14(b). As a function of t, the number rj moves clockwise around 
the circle of radius a with angular velocity co. 

It is important that you get a good feel for the role of the complex constant A = ae lS 
in (2.77): If A happened to equal 1, then rj would be just r] = e~ lC0t , which lies on the 
unit circle, moving clockwise with angular velocity co and starting from the real axis 
{ji = 1) when t — 0. If A = a is real but not equal to 1, then it simply magnifies the 
unit circle to a circle of radius a, around which t] moves with the same angular speed 
and starting from the real axis, at rj = a when t = 0. Finally if A = ae lh , then the 
effect of the angle 8 is to rotate rj through the fixed angle <5, so that ij starts out at t = 0 
with polar angle <5. 

Armed with these mathematical results, we can now return to the charged particle 
in a magnetic field. 
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2.7 Solution for the Charge in a B Field 


Mathematically, the solution for the velocity v of our charged particle in a B field is 
complete, and all that remains is to interpret it physically. We already know that v z , 
the component along B, is constant. The components (v x , v y ) transverse to B we have 
represented by the complex number rj = v x + iv y , and we have seen that Newton’s 
second law implies that q has the time dependence rj = Ae~ im , moving uniformly 
around the circle of Figure 2.14(b). Now, the arrow shown in that figure, pointing 
from O to rj, is in fact a pictorial representation of the transverse velocity (v x , v y ). 
Therefore this transverse velocity changes direction, turning clockwise, with constant 
angular velocity 7 co = qB/m and with constant magnitude. Because v z is constant, 
this suggests that the particle undergoes a spiralling, or helical, motion. To verify this, 
we have only to integrate v to find r as a function of t. 

That v z is constant implies that 

z(t ) = z 0 + v zo t. (2.78) 

The motion of x and y is most easily found by introducing another complex number 

£ = x + iy 


where £ is the Greek letter xi. In the complex plane, c, is a picture of the transverse 
position (x, y). Clearly, the derivative of £ is rj, that is, ^ = rj. Therefore, 



= — e l(ot + constant. (2.79) 

co 

If we rename the coefficient iA/co as C and the constant of integration as X + iY, 
this implies that 

x + iy = Ce~ i0Jt + (X + iY). 

By redefining our origin so that the z axis goes through the point ( X, Y), we can 
eliminate the constant term on the right to give 

x + iy = Ce~ i(ot , (2.80) 

and, by setting t = 0, we can identify the remaining constant C as 

C = x 0 + iy 0 . 

This result is illustrated in Figure 2.15. We see there that the transverse position (x, y) 
moves clockwise round a circle with angular velocity co = qB/m. Meanwhile z as 
given by (2.78) increases steadily, so the particle actually describes a uniform helix 
whose axis is parallel to the magnetic field. 


7 1 am assuming the charge q is positive; if q is negative, then co = qB/m is negative, meaning 
that the transverse velocity rotates counterclockwise. 
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Figure 2.15 Motion of a charge in a uniform magnetic 
field in the z direction. The transverse position (x, y) 
moves around a circle as shown, while the coordinate z 
moves with constant velocity into or out of the page. 


There are many examples of the helical motion of a charged particle along a 
magnetic field; for example, cosmic-ray particles (charged particles hitting the earth 
from space) can get caught by the earth’s magnetic field and spiral north or south 
along the field lines. If the z component of the velocity happens to be zero, then 
the spiral reduces to a circle. In the cyclotron, a device for accelerating charged 
particles to high energies, the particles are trapped in circular orbits in this way. They 
are slowly accelerated by the judiciously timed application of an electric field. The 
angular frequency of the orbit is, of course, co = qB/m (which is why this is called 
the cyclotron frequency). The radius of the orbit is 


v mv p 
co qB qB 


(2.81) 


This radius increases as the particles accelerate, so that they eventually emerge at the 
outer edge of the circular magnets that produce the magnetic field. 

The same method that we have used here for a charge in a magnetic field can 
also be used for a particle in magnetic and electric fields, but since this complication 
adds nothing to the method of solution, I shall leave you to try it for yourself in 
Problems 2.53 and 2.55. 


Principal Definitions and Equations of Chapter 2 _ 

Linear and Quadratic Drags 

Provided the speed v is well below that of sound, the magnitude of the drag force 
f = —f(v)\ on an object moving through a fluid is usually well approximated as 


f(v) = flin + /quad 


where 


f\m = bv = pDv and / quad = cv 2 = yD 2 v 2 . [Eqs. (2.2) to (2.6)] 
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Here D denotes the linear size of the object. For a sphere, D is the diameter and, for 
a sphere in air at STP, fi — 1.6 x 10 “ 4 N-s/m 2 and y — 0.25 N-s 2 /m 4 . 

The Lorentz Force on a Charged Particle 

F = ^(E + vx B). [Eq. (2.62) & Problem 2.53] 


Problems for Chapter 2 _ 

Stars indicate the approximate level of difficulty, from easiest (*) to most difficult (*+*). 


section 2.1 Air Resistance 

2.1 * When a baseball flies through the air, the ratio / quac i//iin °f the quadratic to the linear drag force 
is given by (2.7). Given that a baseball has diameter 7 cm, find the approximate speed v at which the 
two drag forces are equally important. For what approximate range of speeds is it safe to treat the drag 
force as purely quadratic? Under normal conditions is it a good approximation to ignore the linear 
term? Answer the same questions for a beach ball of diameter 70 cm. 

2.2 * The origin of the linear drag force on a sphere in a fluid is the viscosity of the fluid. According 
to Stokes’s law, the viscous drag on a sphere is 


fun = InrjDv 


(2.82) 


where r] is the viscosity 8 9 of the fluid, D the sphere’s diameter, and v its speed. Show that this expression 
reproduces the form (2.3) for / lin , with b given by (2.4) as b — fiD. Given that the viscosity of air at 
STP is r) = 1.7 x 10~ 5 N-s/m 2 , verify the value of given in (2.5). 

2.3 ★ (a) The quadratic and linear drag forces on a moving sphere in a fluid are given by (2.84) and 
(2.82) (Problems 2.4 and 2.2). Show that the ratio of these two kinds of drag force can be written as 
/quad//lin = R/ 48, 9 where the dimensionless Reynolds number R is 


Dvq 

V 


(2.83) 


where D is the sphere’s diameter, v its speed, and q and rj are the fluid’s density and viscosity. Clearly 
the Reynolds number is a measure of the relative importance of the two kinds of drag. 10 When R is 


8 For the record, the viscosity r] of a fluid is defined as follows: Imagine a wide channel along which fluid is 
flowing (.x direction) such that the velocity v is zero at the bottom (y = 0) and increases toward the top (y = /z), so 
that successive layers of fluid slide across one another with a velocity gradient dv/dy . The force F with which an 
area A of any one layer drags the fluid above it is proportional to A and to dv/dy , and r] is defined as the constant 
of proportionality; that is, F = r] A dv/dy. 

9 The numerical factor 48 is for a sphere. A similar result holds for other bodies, but the numerical factor is 
different for different shapes. 

10 The Reynolds number is usually defined by (2.83) for flow involving any object, with D defined as a typical 
linear dimension. One sometimes hears the claim that R is the ratio / qua d//iin- Since / qua d//iin = R /48 for a sphere, 
this claim would be better phrased as “R is roughly of the order of / qua d//iin •” 
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very large, the quadratic drag is dominant and the linear can be neglected; vice versa when R is very 
small, (b) Find the Reynolds number for a steel ball bearing (diameter 2 mm) moving at 5 cm/s through 
glycerin (density 1.3 g/cm 3 and viscosity 12 N s/m 2 at STP). 

2.4 ** The origin of the quadratic drag force on any projectile in a fluid is the inertia of the fluid that the 
projectile sweeps up. (a) Assuming the projectile has a cross-sectional area A (normal to its velocity) 
and speed v, and that the density of the fluid is q, show that the rate at which the projectile encounters 
fluid (mass/time) is qAv. (b) Making the simplifying assumption that all of this fluid is accelerated 
to the speed v of the projectile, show that the net drag force on the projectile is qAv 1 . It is certainly 
not true that all the fluid that the projectile encounters is accelerated to the full speed v, but one might 
guess that the actual force would have the form 

/quad = KQAV 2 (2.84) 

where k is a number less than 1, which would depend on the shape of the projectile, with k small for 
a streamlined body, and larger for a body with a flat front end. This proves to be true, and for a sphere 
the factor k is found to be k = 1/4. (c) Show that (2.84) reproduces the form (2.3) for / quad , with c 
given by (2.4) as c = yD 2 . Given that the density of air at STP is q = 1.29 kg/m 3 and that k = 1/4 for 
a sphere, verify the value of y given in (2.6). 

section 2.2 Linear Air Resistance 

2.5 * Suppose that a projectile which is subject to a linear resistive force is thrown vertically down 
with a speed v yo which is greater than the terminal speed u ter . Describe and explain how the velocity 
varies with time, and make a plot of v y against t for the case that v y0 = 2v ter . 

2.6 * (a) Equation (2.33) gives the velocity of an object dropped from rest. At first, when v y is small, 
air resistance should be unimportant and (2.33) should agree with the elementary result v y = gt 
for free fall in a vacuum. Prove that this is the case. [Hint: Remember the Taylor series for e x = 
1 + x + A' 2 / 2 ! + x 3 /3! + • • •, for which the first two or three terms are certainly a good approximation 
when x is small.] (b) The position of the dropped object is given by (2.35) with v y0 = 0. Show similarly 
that this reduces to the familiar y = \gt 2 when t is small. 

2.7 ★ There are certain simple one-dimensional problems where the equation of motion (Newton’s 
second law) can always be solved, or at least reduced to the problem of doing an integral. One of these 
(which we have met a couple of times in this chapter) is the motion of a one-dimensional particle subject 
to a force that depends only on the velocity v, that is, F = F(v). Write down Newton’s second law and 
separate the variables by rewriting it as m dv/F(v ) = dt. Now integrate both sides of this equation and 
show that 


t = m 



dv' 
F(v ')' 


Provided you can do the integral, this gives f as a function of v. You can then solve to give v as a 
function of t. Use this method to solve the special case that F(v) = F 0 , a constant, and comment on 
your result. This method of separation of variables is used again in Problems 2.8 and 2.9. 

2.8 * A mass m has velocity v 0 at time t = 0 and coasts along the x axis in a medium where the drag 
force is F(v) = — cv 3 ? 2 . Use the method of Problem 2.7 to find v in terms of the time t and the other 
given parameters. At what time (if any) will it come to rest? 
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2.9* We solved the differential equation (2.29), mv y = —b(v y — u ter ), for the velocity of an object 
falling through air, by inspection — a most respectable way of solving differential equations. Never¬ 
theless, one would sometimes like a more systematic method, and here is one. Rewrite the equation in 
the “separated” form 


m dv v 

-= -bdt 

Vy ~ V teT 

and integrate both sides from time 0 to t to find v y as a function of t. Compare with (2.30). 

2.10 ** For a steel ball bearing (diameter 2 mm and density 7.8 g/cm 3 ) dropped in glycerin (density 
1.3 g/cm 3 and viscosity 12 N-s/m 2 at STP), the dominant drag force is the linear drag given by (2.82) 
of Problem 2.2. (a) Find the characteristic time r and the terminal speed v ter . [In finding the latter, 
you should include the buoyant force of Archimedes. This just adds a third force on the right side of 
Equation (2.25).] How long after it is dropped from rest will the ball bearing have reached 95% of its 
terminal speed? (b) Use (2.82) and (2.84) (with /c = 1/4 since the ball bearing is a sphere) to compute 
the ratio / qua d//ii n at the terminal speed. Was it a good approximation to neglect / quad ? 

2.11 ** Consider an object that is thrown vertically up with initial speed u 0 in a linear medium, 
(a) Measuring y upward from the point of release, write expressions for the object’s velocity v y (t) 
and position y(t). (b) Find the time for the object to reach its highest point and its position y max at that 
point, (c) Show that as the drag coefficient approaches zero, your last answer reduces to the well-known 
result y max = \v 2 jg for an object in the vacuum. [Hint: If the drag force is very small, the terminal 
speed is very big, so u 0 /u ter is very small. Use the Taylor series for the log function to approximate 
ln(l + 8) by <5 — )<S 2 . (For a little more on Taylor series see Problem 2.18.)] 

2.12 ** Problem 2.7 is about a class of one-dimensional problems that can always be reduced to doing 
an integral. Here is another. Show that if the net force on a one-dimensional particle depends only on 
position, F = Fix), then Newton’s second law can be solved to find v as a function of x given by 

v 2 = v] + — [ F(x') dx'. (2.85) 

m Jxo 

[Hint: Use the chain rule to prove the following handy relation, which we could call the “v d v /dx rule”: 
If you regard v as a function of x, then 


dv 1 d v 2 
V = v — =-. 

dx 2 dx 


( 2 . 86 ) 


Use this to rewrite Newton’s second law in the separated form md{v 2 ) = 2F(x) dx and then 
integrate from x 0 to x.] Comment on your result for the case that Fix) is actually a constant. (You 
may recognise your solution as a statement about kinetic energy and work, both of which we shall be 
discussing in Chapter 4.) 

2.13 ** Consider a mass m constrained to move on the x axis and subject to a net force F = —kx where 
k is a positive constant. The mass is released from rest at x = x 0 at time t — 0. Use the result (2.85) 
in Problem 2.12 to find the mass’s speed as a function of x; that is, dx/dt = g(x) for some function 
g(x). Separate this as dx/g(x) = dt and integrate from time 0 to t to find x as a function of t. (You 
may recognize this as one way — not the easiest — to solve the simple harmonic oscillator.) 

2.14 *** Use the method of Problem 2.7 to solve the following: A mass m is constrained to move along 
the x axis subject to a force F(v) = —F 0 e v ^ v , where F 0 and V are constants, (a) Find v(t) if the initial 
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velocity is v 0 > 0 at time t = 0. (b) At what time does it come instantaneously to rest? (c) By integrating 
v(t), you can find x(t). Do this and find how far the mass travels before coming instantaneously to rest. 

section 2.3 Trajectory and Range in a Linear Medium 

2.15 * Consider a projectile launched with velocity (v xo , v yo ) from horizontal ground (with x measured 
horizontally and y vertically up). Assuming no air resistance, find how long the projectile is in the air 
and show that the distance it travels before landing (the horizontal range) is 2v xo v yo /g. 

2.16 * A golfer hits his ball with speed v 0 at an angle 0 above the horizontal ground. Assuming that 
the angle 6 is fixed and that air resistance can be neglected, what is the minimum speed u 0 (min) for 
which the ball will clear a wall of height h, a distance d away? Your solution should get into trouble if 
the angle 6 is such that tan 0 < h /d. Explain. What is w 0 (min) if 0 = 25°, d = 50 m, and h = 2 m? 

2.17 ★ The two equations (2.36) give a projectile’s position ( x , y) as a function of t. Eliminate t to 
give y as a function of x. Verify Equation (2.37). 

2.18 * Taylor’s theorem states that, for any reasonable function fix), the value of / at a point (x + 8) 
can be expressed as an infinite series involving / and its derivatives at the point x: 

f(x + 8) = fix ) + f'ix)8 + I f"ix)8 2 + 1 f"'ix)8 3 + • • • (2.87) 

where the primes denote successive derivatives of fix). (Depending on the function this series may 
converge for any increment 5 or only for values of 8 less than some nonzero “radius of convergence.”) 
This theorem is enormously useful, especially for small values of 8, when the first one or two terms of 
the series are often an excellent approximation. 11 (a) Find the Taylor series for ln(l + 3). (b) Do the 
same for cos 8. (c) Likewise sin <5. (d) And e s . 

2.19 * Consider the projectile of Section 2.3. (a) Assuming there is no air resistance, write down the 
position (x, y) as a function of t, and eliminate t to give the trajectory y as a function of x. (b) The 
correct trajectory, including a linear drag force, is given by (2.37). Show that this reduces to your 
answer for part (a) when air resistance is switched off (r and u ter = gx both approach infinity). [Hint: 
Remember the Taylor series (2.40) for ln(l — e).] 

2.20 ** [Computer] Use suitable graph-plotting software to plot graphs of the trajectory (2.36) of a 
projectile thrown at 45°above the horizontal and subject to linear air resistance for four different values 
of the drag coefficient, ranging from a significant amount of drag down to no drag at all. Put all four 
trajectories on the same plot. [Hint: In the absence of any given numbers, you may as well choose 
convenient values. For example, why not take v x0 = v yo = 1 and g = 1. (This amounts to choosing your 
units of length and time so that these parameters have the value 1.) With these choices, the strength 
of the drag is given by the one parameter v ter = r, and you might choose to plot the trajectories for 
u ter = 0.3, 1, 3, and oo (that is, no drag at all), and for times from t = 0 to 3. For the case that u ter = oo, 
you’ll probably want to write out the trajectory separately.] 

2.21 A gun can fire shells in any direction with the same speed v Q . Ignoring air resistance and 
using cylindrical polar coordinates with the gun at the origin and z measured vertically up, show that 


1 1 

For more details on Taylor’s series see, for example, Mary Boas, Mathematical Methods in the Physical Sci¬ 
ences (Wiley, 1983), p. 22 or Donald McQuarrie, Mathematical Methods for Scientists and Engineers (University 
Science Books, 2003), p. 94. 
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the gun can hit any object inside the surface 


z = 



g 

2v 2 

O 



Describe this surface and comment on its dimensions. 


2.22 ★** [Computer] The equation (2.39) for the range of a projectile in a linear medium cannot be 
solved analytically in terms of elementary functions. If you put in numbers for the several parameters, 
then it can be solved numerically using any of several software packages such as Mathematica, Maple, 
and MatLab. To practice this, do the following: Consider a projectile launched at angle 9 above the 
horizontal ground with initial speed v 0 in a linear medium. Choose units such that v 0 = 1 and g = 1. 
Suppose also that the terminal speed u ter = 1. (With v 0 — t> ter , air resistance should be fairly important.) 
We know that in a vacuum, the maximum range occurs at 6 = tz/A « 0.75. (a) What is the maximum 
range in a vacuum? (b) Now solve (2.39) for the range in the given medium at the same angle 9 = 0.75. 
(c) Once you have your calculation working, repeat it for some selection of values of 9 within which 
the maximum range probably lies. (You could try 9 = 0.4,0.5, • • •, 0.8.) (d) Based on these results, 
choose a smaller interval for 9 where you’re sure the maximum lies and repeat the process. Repeat it 
again if necessary until you know the maximum range and the corresponding angle to two significant 
figures. Compare with the vacuum values. 


section 2.4 Quadratic Air Resistance 

2.23 * Find the terminal speeds in air of (a) a steel ball bearing of diameter 3 mm, (b) a 16-pound steel 
shot, and (c) a 200-pound parachutist in free fall in the fetal position. In all three cases, you can safely 
assume the drag force is purely quadratic. The density of steel is about 8 g/cm 3 and you can treat the 
parachutist as a sphere of density 1 g/cm 3 . 

2.24 * Consider a sphere (diameter D, density £ sph ) falling through air (density £ air ) and assume that 
the drag force is purely quadratic, (a) Use Equation (2.84) from Problem 2.4 (with tc = 1/4 for a sphere) 
to show that the terminal speed is 


v 


ter 



Qsph 

£?air 


( 2 . 88 ) 


(b) Use this result to show that of two spheres of the same size, the denser one will eventually fall 
faster, (c) For two spheres of the same material, show that the larger will eventually fall faster. 

2.25 ★ Consider the cyclist of Section 2.4, coasting to a halt under the influence of a quadratic drag 

force. Derive in detail the results (2.49) and (2.51) for her velocity and position, and verify that the 

> 

constant r = m/cv 0 is indeed a time. 

2.26 * A typical value for the coefficient of quadratic air resistance on a cyclist is around c = 0.20 
N/(m/s) 2 . Assuming that the total mass (cyclist plus cycle) is m = 80 kg and that at t = 0 the cyclist 
has an initial speed v 0 = 20 m/s (about 45 mi/h) and starts to coast to a stop under the influence of air 
resistance, find the characteristic time r = m/cv 0 . How long will it take him to slow to 15 m/s? What 
about 10 m/s? And 5 m/s? (Below about 5 m/s, it is certainly not reasonable to ignore friction, so there 
is no point pursuing this calculation to lower speeds.) 
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2.27 * I kick a puck of mass m up an incline (angle of slope = 9) with initial speed v Q . There is no 
friction between the puck and the incline, but there is air resistance with magnitude f(v) = cv 2 . Write 
down and solve Newton’s second law for the puck’s velocity as a function of t on the upward journey. 
How long does the upward journey last? 

2.28 * A mass m has speed v 0 at the origin and coasts along the x axis in a medium where the drag 
force is F(v) = —cv 3 / 2 . Use the “u dv/dx rule” (2.86) in Problem 2.12 to write the equation of motion 
in the separated form m vdv/F(v ) = dx, and then integrate both sides to give x in terms of v (or vice 
versa). Show that it will eventually travel a distance 

2.29* The terminal speed of a 70-kg sky diver in spread-eagle position is around 50 m/s (about 115 
mi/h). Find his speed at times t = 1, 5, 10, 20, 30 seconds after he jumps from a stationary balloon. 
Compare with the corresponding speeds if there were no air resistance. 

2.30 ★ Suppose we wish to approximate the skydiver of Problem 2.29 as a sphere (not a very promising 
approximation, but nevertheless the kind of approximation physicists sometimes like to make). Given 
the mass and terminal speed, what should we use for the diameter of the sphere? Does your answer 
seem reasonable? 

2.31 ** A basketball has mass m = 600 g and diameter D = 24 cm. (a) What is its terminal speed? 

(b) If it is dropped from a 30-m tower, how long does it take to hit the ground and how fast is it going 
when it does so? Compare with the corresponding numbers in a vacuum. 

2.32 ** Consider the following statement: If at all times during a projectile’s flight its speed is much 
less than the terminal speed, the effects of air resistance are usually very small, (a) Without reference 
to the explicit equations for the magnitude of u ter , explain clearly why this is so. (b) By examining the 
explicit formulas (2.26) and (2.53) explain why the statement above is even more useful for the case 
of quadratic drag than for the linear case. [Hint: Express the ratio // mg of the drag to the weight in 
terms of the ratio n/n ter .] 

2.33 ** The hyperbolic functions cosh z and sinh z are defined as follows: 

€ Z + C~ Z C Z — 6 Z 

cosh z =- and sinh z =- 

2 2 

for any z, real or complex, (a) Sketch the behavior of both functions over a suitable range of real 
values of z. (b) Show that cosh z = cos(iz). What is the corresponding relation for sinh z? (c) What are 
the derivatives of cosh z and sinh z? What about their integrals? (d) Show that cosh 2 z — sinh 2 z = 1. 
(e) Show that f dx/y/l + x 2 = arcsinhx. [Hint: One way to do this is to make the substitution 
x = sinhz.] 

2.34** The hyperbolic function tanhz is defined as tanhz = sinhz/coshz, with coshz and sinhz 
defined as in Problem 2.33. (a) Prove that tanhz = —i tan(iz). (b) What is the derivative of tanhz? 

(c) Show that f dz tanhz = In coshz. (d) Prove that 1 — tanh 2 z = sech 2 z, where sechz = 1/ coshz. 
(e) Show that f dx/( 1 — x 2 ) = arctanhx. 

2.35 ★* (a) Fill in the details of the arguments leading from the equation of motion (2.52) to Equations 
(2.57) and (2.58) for the velocity and position of a dropped object subject to quadratic air resistance. 
Be sure to do the two integrals involved. (The results of Problem 2.34 will help.) (b) Tidy the two 
equations by introducing the parameter r = v ter /g. Show that when t = r, v has reached 76% of its 
terminal value. What are the corresponding percentages when t = 2r and 3r? (c) Show that when 
t t, the position is approximately y ~ v ter t + const. [Hint: The definition of coshx (Problem 2.33) 
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gives you a simple approximation when x is large.] (d) Show that for t small, Equation (2.58) for the 
position gives y ~ \gt 2 . [Use the Taylor series for cosh x and for ln(l + 5).] 

2.36 ** Consider the following quote from Galileo’s Dialogues Concerning Two New Sciences: 

Aristotle says that “an iron ball of 100 pounds falling from a height of one hundred cubits reaches the 
ground before a one-pound ball has fallen a single cubit.” I say that they arrive at the same time. You 
find, on making the experiment, that the larger outstrips the smaller by two finger-breadths, that is, when 
the larger has reached the ground, the other is short of it by two finger-breadths. 

We know that the statement attributed to Aristotle is totally wrong, but just how close is Galileo’s claim 
that the difference is just “two finger breadths”? (a) Given that the density of iron is about 8 g/cm, find 
the terminal speeds of the two iron balls, (b) Given that a cubit is about 2 feet, use Equation (2.58) to 
find the time for the heavier ball to land and then the position of the lighter ball at that time. How far 
apart are they? 

2.37 ** The result (2.57) for the velocity of a falling object was found by integrating Equation (2.55) 
and the quickest way to do this is to use the integral f du/{\ — ir ) = arctanh u. Here is another way 
to do it: Integrate (2.55) using the method of “partial fractions,” writing 

1 — u 2 2 \ 1 + w 1 — u ) 

which lets you do the integral in terms of natural logs. Solve the resulting equation to give v as a 
function of t and show that your answer agrees with (2.57). 

2.38** A projectile that is subject to quadratic air resistance is thrown vertically up with initial 
speed v Q . (a) Write down the equation of motion for the upward motion and solve it to give rasa 
function of t. (b) Show that the time to reach the top of the trajectory is 

^top IT ter/ s') arctan(i> 0 /i> ter ). 

(c) For the baseball of Example 2.5 (with u ter = 35 m/s), find f top for the cases that v 0 = 1, 10, 20, 30, 
and 40 m/s, and compare with the corresponding values in a vacuum. 

2.39 ★★ When a cyclist coasts to a stop, he is actually subject to two forces, the quadratic force of air 
resistance, / = — cv 2 (with c as given in Problem 2.26), and a constant frictional force / fr of about 3 
N. The former is dominant at high and medium speeds, the latter at low speed. (The frictional force is a 
combination of ordinary friction in the bearings and rolling friction of the tires on the road.) (a) Write 
down the equation of motion while the cyclist is coasting to a stop. Solve it by separating variables to 
give t as a function of v. (b) Using the numbers of Problem 2.26 (and the value / fr = 3 N given above) 
find how long it takes the cyclist to slow from his initial 20 m/s to 15 m/s. How long to slow to 10 and 
5 m/s? How long to come to a full stop? If you did Problem 2.26, compare with the answers you got 
there ignoring friction entirely. 

2.40 ** Consider an object that is coasting horizontally (positive x direction) subject to a drag force 
/ = —bv — cv 2 . Write down Newton’s second law for this object and solve for v by separating 
variables. Sketch the behavior of v as a function of t. Explain the time dependence for t large. (Which 
force term is dominant when t is large?) 

2.41 ** A baseball is thrown vertically up with speed v 0 and is subject to a quadratic drag with 
magnitude f(v ) = cv 2 . Write down the equation of motion for the upward journey (measuring y 
vertically up) and show that it can be rewritten as v = —g[l + (u/u ter ) 2 ]- Use the “vdv/dx rule” 
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(2.86) to write i) as vdv/dy, and then solve the equation of motion by separating variables (put all 
terms involving v on one side and all terms involving y on the other). Integrate both sides to give y in 
terms of v, and hence v as a function of y. Show that the baseball’s maximum height is 


v 


Tmax 


ter 


2g 


In 


/ v l + v 

\ u 2 

\ ter 



(2.89) 


If v 0 = 20 m/s (about 45 mph) and the baseball has the parameters given in Example 2.5 (page 61), 
what is y max ? Compare with the value in a vacuum. 

2.42** Consider again the baseball of Problem 2.41 and write down the equation of motion for the 
downward journey. (Notice that with a quadratic drag the downward equation is different from the 
upward one, and has to be treated separately.) Find v as a function of y and, given that the downward 
journey starts at y max as given in (2.89), show that the speed when the ball returns to the ground is 

^ter^o/y^ter + v o' Discuss this result for the cases of very much and very little air resistance. What 
is the numerical value of this speed for the baseball of Problem 2.41? Compare with the value in a 
vacuum. 


2.43 *** [Computer] The basketball of Problem 2.31 is thrown from a height of 2 m with initial velocity 
v 0 = 15 m/s at 45° above the horizontal, (a) Use appropriate software to solve the equations of motion 
(2.61) for the ball’s position (x, y) and plot the trajectory. Show the corresponding trajectory in the 
absence of air resistance, (b) Use your plot to find how far the ball travels in the horizontal direction 
before it hits the floor. Compare with the corresponding range in a vacuum. 

2.44 *★* [Computer] To get an accurate trajectory for a projectile one must often take account of several 
complications. For example, if a projectile goes very high then we have to allow for the reduction 
in air resistance as atmospheric density decreases. To illustrate this, consider an iron cannonball 
(diameter 15 cm, density 7.8 g/cm 3 ) that is fired with initial velocity 300 m/s at 50 degrees above 
the horizontal. The drag force is approximately quadratic, but since the drag is proportional to the 
atmospheric density and the density falls off exponentially with height, the drag force is / = c{y)v 2 
where c(y) = yD 2 exp(— y/X) with y given by (2.6) and X ~ 10, 000 m. (a) Write down the equations 
of motion for the cannonball and use appropriate software to solve numerically for x{t) and y(t) for 
0 < t < 3.5 s. Plot the ball’s trajectory and find its horizontal range, (b) Do the same calculation ignoring 
the variation of atmospheric density [that is, setting c(y) = c(0)], and yet again ignoring air resistance 
entirely. Plot all three trajectories for 0 < t < 3.5 s on the same graph. You will find that in this case 
air resistance makes a huge difference and that the variation of air resistance makes a small, but not 
negligible, difference. 


section 2.6 Complex Exponentials 

2.45 * (a) Using Euler’s relation (2.76), prove that any complex number z — x + iy can be written in 
the form z = re’ e , where r and 6 are real. Describe the significance of r and 0 with reference to the 
complex plane, (b) Write z = 3 + 4i in the form z = re'" . (c) Write z = 2e~ in ^ in the form x + iy. 

2.46 * For any complex number z = x + iy, the real and imaginary parts are defined as the real 
numbers Re(z) = x and Im(z) = y. The modulus or absolute value is |zj — yx 2 + y 2 and the phase 
or angle is the value of 9 when z is expressed as z = re ld . The complex conjugate is z* = x — iy. 
(This last is the notation used by most physicists; most mathematicians use z.) For each of the following 
complex numbers, find the real and imaginary parts, the modulus and phase, and the complex conjugate, 
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and sketch z and z* in the complex plane: 

(a) z = 1 + i (b) z = 1 — i V3 
(c) z = «j2e~ in l A (d) z = 5e ia)t . 

In part (d), co is a constant and t is the time. 

2.47 * For each of the following two pairs of numbers compute z + w, z — w, zw, and z/w. 

(a) z = 6 + 8/ and w = 3 — 4i (b) z = Se ,Jt / 3 and w = 4e l7r / 6 . 

Notice that for adding and subtracting complex numbers, the form x + iy is more convenient, but for 
multiplying and especially dividing, the form re 10 is more convenient. In part (a), a clever trick for 
finding z/w without converting to the form re 10 is to multiply top and bottom by in*; try this one both 
ways. 

2.48 * Prove that |z| = \fzFz for any complex number z. 

2.49 * Consider the complex number z = e 10 = cos 9 + i sin 0. (a) By evaluating z 2 two different ways, 
prove the trig identities cos 26 — cos 2 6 — sin 2 6 and sin 26 = 2 sin 6 cos 6. (b) Use the same technique 
to find corresponding identities for cos 3 6 and sin 36. 

2.50 * Use the series definition (2.72) of e z to prove that 12 de z /dz = e z . 

2.51 ★★ Use the series definition (2.72) of e z to prove that e z e w = e z+w . [Hint: If you write down the 
left side as a product of two series, you will have a huge sum of terms like z n w m . If you group together 
all the terms for which n + m is the same (call it p) and use the binomial theorem, you will find you 
have the series for the right side.] 

section 2.7 Solution for the Charge in a B Field 

2.52* The transverse velocity of the particle in Sections 2.5 and 2.7 is contained in (2.77), since 
r? = v x + iv y . By taking the real and imaginary parts, find expressions for v x and v y separately. Based 
on these expressions describe the time dependence of the transverse velocity. 

2.53 * A charged particle of mass m and positive charge q moves in uniform electric and magnetic 
fields, E and B, both pointing in the z direction. The net force on the particle is F = q(E + v x B). 
Write down the equation of motion for the particle and resolve it into its three components. Solve the 
equations and describe the particle’s motion. 

2.54 ** In Section 2.5 we solved the equations of motion (2.68) for the transverse velocity of a charge 
in a magnetic field by the trick of using the complex number q = v x + iv y . As you might imagine, 
the equations can certainly be solved without this trick. Here is one way: (a) Differentiate the first of 
equations (2.68) with respect to t and use the second to give you a second-order differential equation 
for v x . This is an equation you should recognize [if not, look at Equation (1.55)] and you can write 
down its general solution. Once you know v x , (2.68) tells you v y . (b) Show that the general solution 
you get here is the same as the general solution contained in (2.77), as disentangled in Problem 2.52. 


12 If you are the type who worries about mathematical niceties, you may be wondering if it is permissible 
to differentiate an infinite series. Fortunately, in the case of a power series (such as this), there is a theorem 
that guarantees the series can be differentiated for any z inside the “radius of convergence.” Since the radius of 
convergence of the series for e z is infinite, we can differentiate it for any value of z. 
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2.55 *** A charged particle of mass m and positive charge q moves in uniform electric and magnetic 
fields, E pointing in the y direction and B in the z direction (an arrangement called “crossed E and B 
fields”)- Suppose the particle is initially at the origin and is given a kick at time t = 0 along the x axis 
with v x = v x0 (positive or negative), (a) Write down the equation of motion for the particle and resolve 
it into its three components. Show that the motion remains in the plane z = 0. (b) Prove that there is 
a unique value of v xo , called the drift speed u dr , for which the particle moves undeflected through the 
fields. (This is the basis of velocity selectors, which select particles traveling at one chosen speed from 
a beam with many different speeds.) (c) Solve the equations of motion to give the particle’s velocity 
as a function of t, for arbitrary values of v xo . [Hint: The equations for (v x , v y ) should look very like 
Equations (2.68) except for an offset of v x by a constant. If you make a change of variables of the 
form u x — v x — u dr and u v = v y , the equations for ( u x , u y ) will have exactly the form (2.68), whose 
general solution you know.] (d) Integrate the velocity to find the position as a function of t and sketch 
the trajectory for various values of v x0 . 



